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ABSTRACT 

The galaxy merger and accretion rates, and their evolution with time, provide important tests for 
models of galaxy formation and evolution. Close pairs of galaxies are the best available means of 
measuring redshift evolution in these quantities. In this study, we introduce two new pair statistics, 
which relate close pairs to the merger and accretion rates. We demonstrate the importance of correcting 
these (and other) pair statistics for selection effects related to sample depth and completeness. In 
particular, we highlight the severe bias that can result from the use of a flux-limited survey. The first 
statistic, denoted Nc, gives the number of companions per galaxy, within a specified range in absolute 
magnitude. Nc is directly related to the galaxy merger rate. The second statistic, called Lc, gives the 
total luminosity in companions, per galaxy. This quantity can be used to investigate the mass accretion 
rate. Both Nc and Lc are related to the galaxy correlation function ^ and luminosity function (f){M) in 
a straightforward manner. Both statistics have been designed with selection effects in mind. We outline 
techniques which account for various selection effects, and demonstrate the success of this approach 
using Monte Carlo simulations. If one assumes that clustering is independent of luminosity (which is 
appropriate for reasonable ranges in luminosity), then these statistics may be applied to flux- limited 
surveys. 

These techniques are applied to a sample of 5426 galaxies in the SSRS2 redshift survey. This is the 
first large, well-defined low-z survey to be used for pair statistics. Using close (5 kpc < < 20 
h-^ kpc) dynamical {Av < 500 km/s) pairs, we find Yc(-21 < Mb < -18) = 0.0226 ± 0.0052 and 
Lc(-21 < Mb < -18) = 0.0216 ±0.0055 x 10^" H'^Lq at z=0.015. These are the first secure estimates of 
low-redshift pair statistics, and they will provide local benchmarks for ongoing and future pair studies. 
If Nc remains fixed with redshift, simple assumptions imply that ~ 6.6% of present day galaxies with 
—21 < Mb < —18 have undergone mergers since z=l. When applied to redshift surveys of more distant 
galaxies, these techniques will yield the first robust estimates of evolution in the galaxy merger and 
accretion rates. 



1. INTRODUCTION 

Studies of galaxy evolution have revealed surprisingly 
recent changes in galaxy populations. Comparisons of 
present day galaxies with those at moderate {z 0.5) 
and high (z ~ 3) redshift have uncovered trends which 
are often dramatic, and may trace galaxies to the time at 
which they were first assembled into recognizable entities. 
These discoveries have shed new light on the formation of 
galaxies, and have provided clues as to the nature of their 
evolution. At z < 1, the picture that is emerging is one in 
which early type galaxies evolve slowly and passively, while 
late type gala xies become mor e numerous with increasing 
redshift (e.g., Lin et al. 1999] ). At higher redshifts, deep 



tions remain regarding the physical processes driving this 
evolution. Mechanisms that have been postulated include 
galaxy-galaxy mergers, luminosity-dependent luminosity 
evolution, and the existence of a new population of galax- 
ies that has faded by the present epoch (see reviews by 



Koo fc Kron 1992 and Ellis 1997). In this study, we will 



surve ys such as the Hubble Deep Field (Williams et al 



1996) indicate an increase in the cosmic star formation 
2 (e.g., Madau, Pozzetti, and Dickinson 



investigate the relative importance of mergers in the evolu- 
tion of field galaxies. Mergers transform the mass function 
of galaxies, marking a progression from small galaxies to 
larger ones. In addition, mergers can completely disrupt 
their constituent galaxies, changing gas-rich spiral galaxies 
into quiescent ellipticals (e.g., Toomre and Toomre 1972). 
During a collision, a merging system may also go through 
a dramatic transition, with the possible onset of triggered 
star formation and/or accretion onto a central black hole 
(see review by Barnes & Hernquist 1992). 

It is clear that mergers do occur, even during the rela- 
tively quiet present epoch. However, the frequency of these 
events, and the distribution of masses involved, has yet to 
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rate out to z 
1998). 

While considerable progress has been made in the obser- 
vational description of galaxy evolution, important ques- 
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be accurately established. This is true at both low and 
high redshift. Furthermore, while a number of attempts 
have been made, a secure measurement of evolution in 
the galaxy merger rate remains elusive, and a comparable 
measure of the accretion rate has yet to be attempted. 

In this study, we introduce a new approach for relating 
dynamically close galaxy pairs to merger and accretion 
rates. These new techniques yield robust measurements 
for disparate samples, thereby allowing meaningful com- 
paris ons of mergers at low and high redshift. In addition. 



(1997), using a flux-limit ed (B < 14.5 ) sample of galaxies 
from the UGC catalog ( Nilson 1973| ). Using pairs with 
projected physical separations of less than 20 kpc, 
they estimated the local pair fraction to be 4.3 ± 0.4%. 
This result was shown to be consistent with the local pair 
fraction estimates of Carlberg et al. (1994) and Yee & 
EUingson (1995), both of whom also used the UGC cat- 
alog. The pair fraction has been measured for samples 
of galaxies at moderate redshift (z ~ 0.5), yieldi ng pub- 



lished estimates ranging f rom approximately 0% (Woods 



samples, and to studies of both major and minor mergers. 
We apply these techniques to a large sample of galaxies 
at low redshift (SSRS2), providing a much needed local 
benchmark for comparison with samples at higher redshift. 
In a forthcoming paper (Patton et al. 2000), we will apply 
these techniques to a large sample of galaxies at moderate 
redshift (CN0C2; 0.1 < z < 0.6), yielding a secure esti- 
mate for the rate of evolution in the galaxy merger and 
accretion rates. 

An overview of earlier pair studies, and a discussion 
of their limitations and shortcomings, are given in the 
next section. The SSRS2 data are described in § Sec- 
tion ^ discusses the connection between close pairs and the 
merger and accretion rates, while § || introduces new statis- 
tics for relating these quantities. Section ^ describes how 
these statistics can be applied to flux-limited surveys in 
a robust manner. A pair classification experiment is pre- 
sented in § 1^ giving empirical justification for our close 
pair criteria. Pair statistics are then computed for the 
SSRS2 survey in § ^, and the implications are discussed in 
§ ^ Conclusions are given in the final section. Through- 
out this paper, we use a Hubble constant of Hq = lOO/i 
km Mpc~^. We assume h=l and qQ=0.1, unless stated 
otherwise. 



u 

these jair statistics can be adapted to a variety of redshift Fahlman, fc Richer 1995 ) to 34% ± 9% ( [Burkey et al 



2. BACKGROUND 

Every estimate of evolution in the merger and/or ac- 
cretion rate begins with the definition of a merger statis- 
tic. Ideally, this statistic should be independent of selec- 
tion effects such as optical contamination due to unrelated 
foreground/background galaxies, redshift incompleteness, 
redshift-dependent changes in minimum luminosity result- 
ing from flux limits, contamination due to non-merging 
systems, /c-corrections, and luminosity evolution. In ad- 
dition, it should be straightforward to relate the statistic 
to the global galaxy population, and to measurements on 
larger scales. The statistic should then be applied to large, 
well-defined samples from low to high redshift, yielding se- 
cure estimates of how the merger and/or accretion rates 
vary with redshift. 

Within the past decade, there have been a number of 



using close pairs of galaxies (e.g., ^epf & Koo 1989|, 


Burke;y 


et al. 1994, Carlberg, Pritchet, & Infante 1994, 


Woods, 


Fahlman, & Richer 1995, Yee & EUingson 199£, 


Pattor 


et al. 1997], Le Fcvrc ct al. 1999). The statistic that 



has been most commonly employed is the traditional pair 
fraction, which gives the fraction of galaxies with suit- 
ably close physical companions. This statistic is assumed 
to be proportional to the galaxy merger rate. The local 
(low-redshift) pair fraction was estimated by Patton et al. 



1994). Evolution in the galaxy merger rate is often pa- 
rameterized as (1 + z)™. Close pair studies have yielded 
a wide variety of results, spanning the range ^ to ^ 5. 
There are several reasons for the large spread in results. 
First, different methods have been used to relate the pair 
fraction to the merger rate. In addition, some estimates 
have been found to suffer from biases due to optical con- 
tamination or redshift completeness. After taking all of 
these effects into account, Patton et al. (1997) demon- 
strated that most results are broadly consistent with their 
estimate of to = 2.8 ± 0.9, made using the largest redshift 
sample (545 galaxies) to date. 

While this convergence seems promising, all of these re- 
sults have suffered from a number of very significant diffi- 
culties. The central (and most serious) problem has been 
the comparison between low and moderate redshift sam- 
ples. Low-z samples have been poorly defined, due to a 
lack of suitable redshift surveys. In addition, the pair frac- 
tion depends on both the clustering and mean density of 
galaxies. The latter is very sensitive to the limiting ab- 
solute magnitude of galaxies, leading to severe redshift- 
dependent biases when using flux-limited galaxy samples. 
These biases have not been taken into account in the com- 
putation of pair fractions, or in the comparison between 
samples at different redshifts. 

While these problems are the most serious, there are 
several other areas of concern. A lack of redshift informa- 
tion has meant dealing with optical contamination due to 
unrelated foreground and background galaxies. Moreover, 
while one can statistically correct for this contamination, it 
is still not possible to discern low velocity companions from 
those that are physically associated but unbound, unless 
additional redshift information is available. Finally, there 
is no direct connection between the pair fraction and the 
galaxy correlation function (CF) and luminosity function 
(LF), making the results more difficult to interpret. 

To address these issues, we have developed a novel ap- 
proach to measuring pair statistics. We will introduce new 
statistics that overcome many of the afflictions of the tra- 
ditional pair fraction. We will then apply these statistics 
to a large, well-defined sample of galaxies at low redshift. 

3. DATA 



The Se cond Southern Sky Redshift Survey (da Costa et 
al. 1998; hereafter SSRS2) consists of 5426 galaxies with 



TTiB < 15.5, in two regions spanning a total of 1.69 steradi- 
ans in the southern celestial hemisphere. The first region, 
denoted SSRS2 South, has boundaries -40° <d< -2.5° 
and bii < -40°. The second region, SSRS2 North, is 
a more recent addition, and is bounded hy 6 < 0° and 
bji > 35°. Galaxies were selected primarily from the list 
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of non-stellar objects in the Hubble Space Telescope Guide 
Star Catalog, with positions accurate to ~ 1" and pho- 
tometry with an rms scatter of ~ 0.3 magnitudes (Alonso 
et al. 1993, Alonso et al. 1994). Steps were taken to 
ensure that single galaxies were not mistakenly identified 
as clo se pairs, du e to the presence of dust lanes, etc. ( |da| 
Costa et al. 1998). In addition, careful attention was paid 



to cases where a very close pair might be mistaken for a 
single galaxy. This was found to make a negligible contri- 
bution to the catalog as a whole (< 0.1% of galaxies are 
affected). The effect on the pairs analysis in this paper is 
further reduced by imposing a minimum pair separation 
of 5 h-^ kpc (see §|7l). 

The sample now includes redshifts for all galaxies 
brighter than niB < 15.5. We correct all velocities to the 
local group barycenter using Equation 6 from Courteau 
and van den Bergh (1999). We restrict our analysis to the 
redshift range 0.005 < z < 0.05. This eliminates nearby 
galaxies, for which recession velocities are dominated by 
peculiar velocities, giving poor distance estimates. We 
also avoid the sparsely sampled high redshift regime. This 
leaves us with a well-defined sample of 4852 galaxies. 

4. THE GALAXY MERGER RATE AND ACCRETION RATE 
4.1. Definitions 

The primary goal of earlier close pair studies has been 
to determine how the galaxy merger rate evolves with red- 
shift. The merger rate affects the mass function of galax- 
ies, and may also be connected to the cosmic star forma- 
tion rate. Before attempting to measure the merger rate, 
it is important to begin with a clear definition of a merger 
and a merger rate. Here, we refer to mergers between two 
galaxies which are both above some minimum mass or lu- 
minosity. If this minimum corresponds roughly to a typical 
bright galaxy (L*), this criterion can be thought of as se- 
lecting so-called major mergers. We consider two merger 
rate definitions. First, it is of interest to determine the 
number of mergers that a typical galaxy will undergo per 
unit time. In this case, the relevant rate may be termed the 
galaxy merger rate (hereafter T^-mg). A related quantity is 
the total number of mergers taking place per unit time per 
unit co-moving volume. We will refer to this as the volume 
merger rate (hereafter T^-mg^,). Clearly, TZmg^ — nolZ^^, 
where rip is the co-moving number density of galaxies. 

While both of these merger rates provide useful mea- 
sures of galaxy interactions, they have their limitations. 
As one probes to faint luminosities, one will find an in- 
creasing number of faint companions; hence, the number 
of inferred mergers will increase in turn. For all realistic 
LFs, this statistic will become dominated by dwarf galax- 
ies. In addition, it is of interest to determine how the mass 
of galaxies will change due to mergers. To address these 
issues, we will also investigate the rate at which mass is 
being accreted onto a typical galaxy. This quantity, the 
total mass accreted per galaxy per unit time, will be re- 
ferred to as the galaxy accretion rate (hereafter 7?.ac). This 
is related to the rate of mass accretion per unit co-moving 
volume (TZucv) by TZucv = noTZa_c- The mass (or luminos- 
ity) dependence of the accretion rate means that it will 
be dominated by relatively massive (or luminous) galax- 
ies, with dwarfs playing a very minor role unless the mass 
function is very steep. 



4.2. Observable Quantities 

In order to determine TZmg observationally, one may be- 
gin by identifying systems which are destined to merge. 
By combining information about the number of these sys- 
tems and the timescale on which they will undergo merg- 
ers, one can estimate an overall merger rate. Specifically, 
if one identifies N„i ongoing mergers per galaxy, and if the 
average merging timescale for these systems is Tmg, then 
Ti-mg — Njn/Tyng- If the mass involved in these mergers 
(per galaxy) is M^, then T^ac = M^/T^g. 

In practice, direct measurement of these quantities is a 
daunting task. It is difficult to determine if a given system 
will merge; furthermore, estimating the merger timescale 
for individual systems is challenging with the limited infor- 
mation generally available. However, if one simply wishes 
to determine how the merger rate is changing with red- 
shift, then the task is more manageable. If one has the 
same definition of a merger in all samples under consid- 
eration, then it is reasonable to assume that the merger 
timescale is the same for these samples. In this case, we 
are left with the task of measuring quantities which are 
directly proportional to the number or mass of mergers 
per galaxy or per unit co-moving volume. If one wishes to 
consider luminosity instead of mass, the relation between 
mass and luminosity must either be the same at all epochs, 
or understood well enough to correct for the differences. 

We have considered several quantities that fit this de- 
scription. All involve the identification of close physical 
associations of galaxies. A "close companion" is defined 
as a neighbour which will merge within a relatively short 
period of time (Tmg ^ iff), which allows an estimate of 
the instantaneous merger/ accretion rate. If a galaxy is 
destined to undergo a merger in the very near future, it 
must have a companion close at hand. One might attempt 
to estimate the number of mergers taking place within a 
sample of galaxies. For example, a close pair of galax- 
ies would be considered one merger, while a close triple 
would lead to two mergers, etc. Owing to the difficulty of 
determining with certainty which systems are undergoing 
mergers, we will not use this approach. One alternative is 
to estimate the number of galaxies with one or more close 
companions, otherwise known as the pair fraction. One 
drawback of this approach is that close triples or higher 
order N-tuples complicate the analysis, since they are re- 
lated to higher orders of the correlation function. This also 
makes it difficult to correct for the flux-limited nature of 
most redshift surveys. As a result, we choose to steer clear 
of this method also. 

In this study, we choose instead to use the number and 
luminosity of close companions per galaxy. The number 
of close companions per galaxy, hereafter Nc, is similar in 
nature to the pair fraction. In fact, they are identical in 
a volume-limited sample with no triples or higher order 
N-tuples. However, Nc will prove to be much more robust 
and versatile. We assume that Nc is directly proportional 
to the number of mergers per galaxy, such that Nm = kNc 
{k is a constant). This pairwise statistic is preferable to 
the number of mergers per galaxy or the fraction of galax- 
ies in merging systems, in that it is related, in a direct and 
straightforward manner, to the galaxy two-point CF and 
the LF (see Section ||). We note that it is not necessary 
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that there be a one-to-one correspondence between com- 
panions and mergers, as long as the correspondence is the 
same, on average, in aU samples under consideration. 

Using this approach to estimate the number of merg- 
ers per galaxy, the merger rate is then given by TZmg — 
kNc/T^g. The actual value of k depends on the merging 
systems under consideration. If one identifies a pure set 
of galaxy pairs, each definitely undergoing a merger, then 
each pair, consisting of 2 companions, would lead to one 
merger, giving A:=0.5. For a pair sample which includes 
some triples and perhaps higher order N-tuples, k < 0.5. 
If the merging sample under investigation contains some 
systems which are not truly merging (for instance, close 
pairs with hyperbolic orbits), then k will also be reduced. 
While k clearly varies with the type of merging system 
used, the key is for k to be the same for all samples under 
consideration. 

We take a similar approach with the accretion rate. We 
again use close companions, and in this case we simply add 
up the luminosity in companions, per galaxy (Lc). Defin- 
ing the mean companion mass-to- light ratio as T, it follows 
that Mm = ^Lc and TZa,c = '^Lc/T^g. When comparing 
different samples, any significant differences in T must be 
accounted for. 

4.3. A Simple Model of Mass Function Evolution Due to 

Mergers 

In order to motivate further the need for merger rate 
measurements, and to set the stage for future work re- 
lating pair statistics to the mass and luminosity function, 
we develop a simple model which relates these important 
quantities. Suppose the galaxy mass function is given by 
4){M,t). This function gives the number density of galax- 
ies of mass M at time t, per unit mass. The model that 
follows can also be expressed in terms of luminosity or 
absolute magnitude, rather than mass. 

We begin by assuming that all changes in the mass func- 
tion are due to mergers. While this is clearly simplistic, 
this model will serve to demonstrate the effects that vari- 
ous merger rates can have on the mass function. In order 
to relate the mass function to the observable luminosity 
function, we further assume that mergers do not induce 
star formation. Again, this is clearly an oversimplifica- 
tion; however, this simple case will still provide a useful 
lower limit on the relative contribution of mergers to LF 
evolution. Finally, we assume that merging is a binary 
process. 



Following Bahcall & Tremaine (1988), we consider how 
4>iM, t) evolves as the universe ages from time t to time 
t+At. Each merger will remove two galaxies from the mass 
function, and produce one new galaxy. Let A0(Af, Ai)out 
represent the decrease in the mass function due to galaxies 
removed by mergers, while A(f>{AI, Ai)in gives the increase 
due to the remnants produced by these mergers. Evolution 
in the mass function can then be given by 

<j>{M, t + At) = 0(Af, t) - A<j>{M, At)out + A(/.(Af, Ai)i„. 

We model this function by considering all galaxy pairs, 
along with an expression for the merging likelihood of each. 
Let p(M, M' , At) denote the probability that a galaxy of 
mass M will merge with a galaxy of mass M' in time in- 
terval At. In order to estimate A(j){M, Ai)out, we need to 



take all galaxies of mass M, and integrate over all com- 
panions, yielding 



A0(M, At)out = 
(j){M, t)(l){M', t)p{M, M', At)[h~^MpcfdM'. (2) 



We devise a comparable expression for A0(Af, At){a by 
integrating over all pairs with end— products of mass M . 
This is achieved by considering all pairs with component 
of mass M - M' and M', such that 



A^(Af, At)i 



(j){M - M',t)(f>{M',t) 



<p{M - M', M', At)[h-^MpcfdM'. 



(3) 



We can also express Equa tion Q in terms of the pair 
statistics outlined in Section |4.2| . If one considers close 
companions of mass M' < oo next to primary galaxies 
of mass M, the volume merger rate can be expressed as 
Unig (M), yielding 



A(/.(Af, At)out = 7^„lg^(Af)A^. 



(4) 



Similarly, if one defines a merger remnant statistic, 
'^mrv(^); to be the co-moving number density of merger 
remnants per unit time corresponding to these same merg- 
ers, then 



A^{M,At)i,,^n^,^{M)At. 



(5) 



Therefore, it is possible, in principle, to use pair statis- 
tics to measure the evolution in the mass or luminosity 
function due to mergers. However, current pair samples 
are too small to permit useful pair statistics for different 
mass combinations. In addition, present day observations 
of close pairs are not of sufficient detail to determine the 
proportion of pairs that will result in mergers (factor k 
in previous section). Moreover, timescale estimates for 
these mergers are not known with any degree of certainty. 
Hence, useful observations of mass function evolution due 
to mergers will have to wait for improved pair samples and 
detailed estimates of merger timescales. 



5. A NEW APPROACH TO MEASURING PAIR STATISTICS 

In this section, we outline the procedure for measuring 
the mean number (Nc) and luminosity (Lc) of close com- 
panions for a sample of galaxies with measured redshifts. 
We begin by defining these statistics in real space, demon- 
strating how they are related to the galaxy LF and CF. We 
then show how these statistics can be applied in redshift 
space. 

5.1. Pair Statistics in Real Space 

In this study, we will measure pair statistics for a com- 
plete low- redshift sample of galaxies (SSRS2). However, 
we wish to make these statistics applicable to a wide vari- 
ety of redshift samples. We would also like this method to 
be useful for studies of minor mergers, where one is inter- 
ested in faint companions around bright galaxies. More- 
over, these techniques should be adaptable to redshift 
samples with varying degrees of completeness (that is, 
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with redshifts not necessarily available for every galaxy). 
Therefore, in the following analysis, we treat host galaxies 
and companions differently 

Consider a primary sample of Ni host galaxies with ab- 
solute magnitudes M < Mi, lying in some volume V. 
Suppose this volume also contains a secondary sample of 
A''2 galaxies with M < M2- In the general case, the pri- 
mary and secondary samples may have galaxies in com- 
mon. This includes the special case in which the two sam- 
ples are identical. If Mi « M2, this will tend to probe 
major mergers. If M2 is chosen to be significantly fainter 
than Ml, this will allow for the study of minor mergers. 
We assume here that both samples are coniplete to the 
given absolute magnitude limits; in Section ra, we extend 
the analysis from volume-limited samples to those that are 
flux-limited. 

We wish to determine the mean number and luminosity 
of companions (in the secondary sample) for galaxies in 
the primary sample. In real space, we define a close com- 
panion to be one that lies at a true physical separation 
of r < r™'*'', where r™'*'^ is some appropriate maximum 
physical separation. To compute the observed mean num- 
ber {Nc) and luminosity (Lc) of companions, we simply 
add up the number {N^ ) and luminosity {La ) of compan- 
ions for each of the Ni galaxies in the primary sample, and 
then compute the mean. Therefore, 



Nc{M < M2) 



and 



E 



NcjNi 



L,(M<M2)=^Le.M. 



(6) 



(7) 



1=1 



We can also estimate what these statistics should be, 
given detailed knowledge of the galaxy two-point CF ^ 
and the LF (j>{M). This is necessary if one wishes to re- 
late these pair statistics to measurements on larger scales. 
Consider a galaxy in the primary sample with absolute 
magnitude Mi at redshift Zi. We would first like to 
estimate the number and luminosity of companions ly- 
ing in a shell at physical distance (proper co-ordinates) 
f l£ fij < 7" + dr from this primary galaxy. To make this 
estimate, we need to know the mean density of galaxies 
(related to the LF), and the expected overdensity in the 
volume of interest (given by the CF) . The mean physical 
number density of galaxies at redshift Zi in the secondary 
sample, with absolute magnitudes M < Mj < M + dM, is 
given by 



n2{zi,M)dM = (1 + Zif(j){M, Zi)dM, 



(8) 



where 4>{M, Zi) is the differential galaxy LF, which spec- 
ifies the co-moving number density of galaxies at redshift 
Zi, in units of /i'^Mpc~'^mag^^. The actual density of ob- 
jects in the region of interest is determined by multiplying 
the mean density by (1-1-^) , where ^ is the overdensity 
given by the two-point CF (Peebles 1980). In general, ^ 
depends on the pair separation r, the mean redshift Zi, 
the absolute magnitude of each galaxy {Mi, Mj), and the 
orbits involved, specified by components parallel (fy) and 
perpendicular {v±_) to the line of sight. It follows that the 



mean number of companions with M < Mj < M -\- dM 
and r < rij < r + dr is given by 

Na{zi,M„M, r, v±,vn)dMdr = 
n2{z„ M)dM[l + e(r, z„ Mi, M, v±,vii)]'iTTr^dr. (9) 

We must now integrate this expression for all compan- 
ions with M < M2 and r < r™^^. Integration over the LF 
yields 

0(M, z,)dM. (10) 

-00 

Integration over the CF is non-trivial, because of the com- 
plex nature of £,{r, Zi, Mi, AI,v±,v\\). With redshift sam- 
ples that are currently available, it is not possible to mea- 
sure this dependence accurately for the systems of inter- 
est. Hence, we must make three important assumptions at 
this stage. First, we assume that ^ is independent of lu- 
minosity. Later in the paper, we demonstrate empirically 
that this is a reasonable assumption, provided one selects a 
sample with appropriate ranges in absolute magnitude (see 
Section 6.1). Secondly, we assume that the distribution of 
velocities is isotropic. If one averages over a reasonable 
number of pairs, this is bound to be true, and therefore 
^ is independent of v±_ and W|| . Finally, we assume that 
the form of the CF, as measured on large scales, can be 
extrapolated to the small scales of interest here. This as- 
sumption applies only to the method of relating pairs to 
large scale measures, and not to the actual measurement 
of pair statistics. 

It is now straightforward to integrate Equation |^. The 
mean number of companions with M < M2 and r < r'^^^ 
for a primary galaxy at redshift Zi is given by 

Nc,{z,,M < M2,r < r'"^'^) = 
n2{z^, M < M2) f [1 + ^{r, z,)]ATTr^dr. (11) 

We derive an analogous expression for the mean luminos- 
ity in companions. The integrated luminosity density is 
given by 

I'M 2 

j2{z„ M < M2) = (1 + z,f / ^{M, Zi)L{M)dM, (12) 

J —OQ 

where 



Therefore, 



L{M) = lO-O-4(*^-*^e)L0. 



LcM.M <M2,r<r'^-^) 



(13) 



j2{z„ M < M2) / [1 + ^(r, z,)]47rr2dr. (14) 

Given measurements of the CF on large scales, it is then 
straightforward to integrate these equations to arrive at 
predicted values of iVc- and Lc^ ■ 

It is important to note that these statistics are directly 
dependent on M2, which affects the mean density of galax- 
ies in the secondary sample. This is different from statis- 
tics such as the correlation function, which are indepen- 
dent of density. Hence, this serves as a reminder that we 
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must exercise caution when choosing our samples, to en- 
sure that differences in the pair statistics (and hence in the 
merger and accretion rates) are not simply due to appar- 
ent density differences resulting from selection effects. In 
addition, note that the choice of A/i has no density-related 
effects on Nr and Lr- 



5.2. Dynamical Pairs in Redshift Space 

While it is preferable to identify companions based on 
their true physical pair separation, this is clearly not fea- 
sible when dealing with data from redshift surveys. In the 
absence of independent distance measurements for each 
galaxy, one must resort to identifying companions in red- 
shift space. In this section, we outline a straightforward 
approach for measuring our new pair statistics in redshift 
space. We then attempt to relate these statistics to their 
counterparts in real space. 

For any given pair of galaxies in redshift space, one can 
compute two basic properties which describe the intrinsic 
pair separation : the projected physical separation (here- 
after Tp) and the rest-frame relative velocity along the line 
of sight (hereafter Av). For a pair of galaxies with redshifts 
Zi (primary galaxy) and Zj (secondary) , with angular sep- 
aration 9, these quantities are given by Vp = 9dA{zi) and 
Av = c\zj — Zi\/{1 + Zi), where dA{zi) is the angular 
diameter distance at redshift Zi. Note that Vp gives the 
projected separation at the redshift of the primary galaxy. 

We define a close companion as one in which the sep- 
aration (both projected and line-of-sight) is less than 
some appropriate separation, such that Vp < r^'^^ and 
Av < Aw™^'^. The line-of-sight criterion depends on both 
the physical line-of-sight separation and the line-of-sight 
peculiar velocity of the companion. It is of course not pos- 
sible to determine the relative contributions of these com- 
ponents without distance information. However, for the 
small companion separations we will be concerned with, 
the peculiar velocity component is likely to be dominant 
in most cases, as we will be dealing with a field sample of 
galaxies (the same would not be true in the high velocity 
environment of rich clusters). Hence, this criterion serves 
primarily to identify companions with low peculiar veloc- 
ities. While this is fundamentally different from the pure 
separation criterion used in real space, it too will serve to 
identify companions with the highest likelihood of under- 
going imminent mergers. Using this definition of a close 
companion, it is straightforward to compute Nc and Lc, 
using Equations ^ and ^. Thus, the complexities of red- 
shift space do not greatly complicate the computation of 
these pair statistics. 

As in real space, we wish to relate these statistics to 
measurements on larger scales, given reasonable assump- 
tions about the LF and CF. The situation is more compli- 
cated in redshift space, and therefore involves additional 
assumptions. We stress, however, that these assumptions 
apply only to the method of relating pair statistics to large 
scale measures, and not to the measured pair statistics 
themselves. To outline an algorithm for generating these 
predictions, we follow the approach of the previous section. 
We begin by modifying Equations |ll| and integrating 
over the new pair volume defined in redshift space. In 
order to do this, we use the two dimensional correlation 



function in redshift space, ^2{rp, r^, z), giving 

n2{zi, M < M2) j J [1 + ^2irp,ry, Zi)]2TTrpdrpdr.v 



and 



L,.(z„M < M2,rp < r-'^^ |r,| < r^^ 



j2 {zi, M < M2) / [1+6 {rp , r„ , Zi)]2Trrpdrpdry . 

The two dimensional correlation function is the convo- 
lution of the velocity distribution in the redshift direction, 
Fivz), with the spatial correlation function ^(r, z), given 

by 



C2(xp,a;„,z) = / S^{Jxl + y'^,z)F{H{z)[y- x^])dy. 

J —OQ 

(17) 

Here, H{z) is the Hubble constant at redshift z, given by 
H{z) = Hq{1 + z)y/T+Tloz. We have ignored the effect 
of infall velocities, which must be taken into account at 
larger radii but is an acceptable approximation for small 
separations. If the form of the CF and LF are known, it is 
straightforward to integrate Equations |l^ and ^ yielding 
predictions of Nc and Lc- 

5.3. Physical Pairs in Redshift Space 

It is not always possible to have precise redshifts for 
all galaxies of interest in a sample. A common scenario 
with redshift surveys is to have redshifts available for a 
subset of galaxies identified in a flux-limited photomet- 
ric sample. The photometric sample used to select galax- 
ies for follow-up spectroscopy probes to fainter apparent 
magnitudes than the spectroscopic sample. In addition, 
the spectroscopic sample may be incomplete, even at the 
bright end of the sample. In this section, we will describe 
the procedure for applying pair statistics to this class of 
samples. 

Suppose the primary sample is defined as all galax- 
ies in the spectroscopic sample with absolute magnitudes 
M < Ml. The secondary sample consists of all galaxies 
lying in the photometric sample, regardless of whether or 
not they have measured redshifts. Once again, there may 
be some overlap between the primary and secondary sam- 
ples. We must now identify close pairs. For each primary- 
secondary pair, we can compute rp in precisely the same 
manner as before (see previous section) , since we need only 
the redshift of the primary galaxy and the angular separa- 
tion of the pair. However, we are no longer able to compute 
the relative velocity along the line of sight, since this re- 
quires redshifts for both members of the pair. Thus, we do 
not have enough information to identify close dynamical 
pairs. However, it is still possible to determine, in a statis- 
tical manner, how many physically associated companions 
are present. This is done by comparing the number (or 
luminosity) of observed companions with the number (or 
luminosity) expected in a ra ndom distribution. 

As stressed in Section 5.1, pair statistics depend on the 



minimum luminosity M2 imposed on the secondary sam- 
ple. While we are now unable to compute the actual lumi- 
nosity for galaxies in the secondary sample, we must still 
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impose M2 if the ensuing pair statistics are to be meaning- 
ful. To do this, we make use of the fact that aU physical 
companions must lie at approximately the same redshift 
as the primary galaxy under consideration. Therefore, M2 
corresponds to a limiting apparent magnitude m2 at red- 
shift Zi, such that 

TO2 ^ M2{zi) - 5\ogh + 25 + 5\ogdL{zi) + k{zi), (18) 

where dL^Zi) is the luminosity distance at redshift z^, and 
k{zi) is the /c-correction. 

To begin, one finds all observed close companions with 
m < 7712 (zi), using only the rp criterion. This results in the 
quantities and L^, where the "D" superscript denote 
companions found in the data sample. One must then es- 
timate the number {N^) and luminosity (L^) of compan- 
ions expected at random. The final pair statistics for close 
physical companions are then given by — — N]^ 
and Lc = L^- Lf. 

We will now describe how to predict these statistics us- 
ing the known LF and CF. This is relatively straightfor- 
ward, since the excess ^ given by the CF is determined by 
the relative proportions of real and random companions. 
The pair statistics are once again integrals over the two 
dimensional CF in redshift space, as specified by Equa- 
tions and ^ In the "1-1- ^" term, the first part gives 
the random contribution, while the second gives the excess 
over random. Thus, these pair statistics give the true den- 
sity of companions, rather than the "excess" density. This 
is intentional, since mergers will occur even in an uncor- 
related, randomly distributed sample of galaxies. At the 
small separations of interest, usually less than 1% of the 
correlation length, the difference between the mean den- 
sity and the mean overdensity is less than about 0.01% in 
real space. For practical measurements in redshift space, 
where r„ = AvH{z)~^ is of order the correlation length, 
the background contribution is substantially larger than 
real space, but still amounts to less than 1%. Thus, for 
the close pairs considered in this study, it is reasonable 
to ignore the contribution that random companions make 
to the sample of physical companions. That is, we take 
1 + ^ ~ This allows us to relate the predictions to 
the measured pair statistics set out above. In principle. 



Equations and 
—00 < r„ < 00 to o 



16| can be integrated over the range 
jtain predictions of and Lc- 



5.4. Application to Volume-limited Monte Carlo 
Simulations 

To illustrate the concepts introduced so far, and to 
emphasize how these statistics depend on M2, we apply 
these techniques to volume-limited Monte Carlo simula- 
tions, which mimic the global distribution of galaxies in 
the SSRS2 North and South catalogs. Using qo—0.5, galax- 
ies are distributed randomly within the co-moving volume 
enclosed by 0.005 < z < 0.05 and the SSRS2 boundaries 
on the sky (see Section ||) . All peculiar velocities are set 
to zero. To create a volume-limited sample, we impose a 
minimum luminosity of Af g=— 16, and assign l uminosi- 
ties using the SSRS2 LF ( [Marzke et al. 1998| ), which 
has Schechter function parameters — 5 log /i^ — 19.43, 
a=-1.12, and = 0.0128 /i^Mpc'^mag"^ An arbitrar- 
ily large number of galaxies can be generated, which is of 



great assistance when looking for small systematic effects. 
We produce 16000 galaxies in the South, and 8070 in the 
North; this gives the same density of galaxies in both re- 
gions. 

Using these simulations, we compute Nc and Lc- As 
these galaxies are distributed randomly (as opposed to 
real galaxies which are clustered), close pairs are relatively 
rare. To ensure a reasonable yield of pairs, we use a pair 
definition of r^>'^'^=l h-^ Mpc and Aw™'^'^=1000 km/s. 
We note that there are no peculiar velocities in these sim- 
ulations; hence, the Av™^^^ criterion provides upper and 
lower limits on the line-of-sight distance to companions. 
Also, recall from the preceding section that the choice of 
Ml has no effect on the pair statistics if clustering is inde- 
pendent of luminosity. Hence, we choose Mi=M2, which 
maximizes the size of the primary sample, and therefore 
minimizes the measurement errors in Nc and Lc- 

With these assumptions, we compute pair statistics for 
a range of choices of Al2. Errors are computed using the 
Jackknife technique. For this resampling method, partial 
standard deviations. Si, are computed for each object by 
taking the difference between the quantity being measur- 
ing, /, and the same quantity with the i*"^ galaxy removed 
from the sample, /i, such that Si = f — fi- For a sample of 
N galaxies, the variance is given by [{N — l)/A^X]i ^iY^^ 
(Efron 1981; Efron & Tibshirani 1986). Results are given 
in Figure ^ Both statistics continue to increase as M2 be- 
comes fainter. Nc diverges at faint magnitudes, while Lc 
is seen to converge. This behaviour is a direct consequence 
of the shape of the LF; Nc converges for a > —1, while 
Lc converges for a > —2. The existence and magnitude of 
these trends clearly demonstrate the need to specify M2 
when measuring pair statistics. 

6. APPLICATION TO FLUX-LIMITED SAMPLES 

The preceding section gives a straightforward prescrip- 
tion for computing pair statistics in volume-limited sam- 
ples. However, redshift surveys are generally flux-limited. 
By defining a volume-limited sample within such a survey, 
one must discard a large proportion of the data. In this 
section, we will outline how these pair statistics can be 
applied to flux-limited surveys. 

Pair statistics necessarily depend on both clustering and 
mean density, as shown by Equations 15 and |l^. In a flux- 
limited sample, both clustering and mean density will vary 
throughout the sample. We will use these equations to ac- 
count for rcdshift-dependent changes in mean density, and 
we will demonstrate how to minimize the effects of cluster- 
ing differences. These techniques will then be tested with 
Monte Carlo simulations. 

6.1. Dependence on Clustering 

By removing the fixed luminosity limit, the overall dis- 
tribution of galaxy luminosities will vary with redshift 
within the sample, and the mean luminosity of the sam- 
ple will differ from the volume-limited sample. However, 
galaxy clustering is known to be luminosity dependent. 
Measures of the galaxy correlation function (e.g., Loveday 
et al. 1995, Willmer et al. 1998), power spectrum (e.g., 
Vogeley 1993), and counts in cells (Benoist et al. 1996) all 
find that luminous galaxies {L > L*) are more clustered 
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than sub-L* galaxies, typically by a factor of ~ 2. This 
increase in clustering may be particularly strong (factor 
^ 4) for very luminous galaxies {Mb < —21). 

Clearly, this effect should not be ignored when comput- 
ing pair statistics. In principle, this could be incorporated 
into the measurement of these pair statistics. However, 
available pair samples are too small to measure this de- 
pendence. We choose instead to minimize these effects by 
restricting the analysis to a fixed range in absolute mag- 
nitude, within which luminosity-dependent clustering is 
small or negligible. This is done by imposing additional 
bright (Mbright) and faint (Mfaint) absolute magnitude lim- 
its on the sample. Having thereby reduced the effects of 
luminosity segregation, we then assume that the remaining 
differences will not have a signifi cant effect on the mea- 
sured pair statistics. In Section 8.3, we use the SSRS2 



sample to demonstrate empirically that this is in fact a 
reasonable assumption. 

6.2. Dependence on Limiting Absolute Magnitude 

In Section ^ we demonstrated that these pair statistics 
are meaningful only if one specifies the minimum lumi- 
nosity of the primary and secondary samples. For a flux- 
limited sample, however, the minimum luminosity of the 
sample increases with redshift. One must therefore decide 
on a representative minimum luminosity, and account for 
differences between the desired minimum luminosity and 
the redshift-dependent minimum imposed by the appar- 
ent magnitude limit of the sample. If the LF is known, 
this can be achieved by weighting each galaxy appropri- 
ately. In this section, we outline a weighting scheme which 
makes this correction. 

6.2.1. Weighting of Secondary Sample 

Consider a flux-limited sample in which host galaxies are 
located at a variety of redshifts. Those at low redshift will 
have the greatest probability of having close companions 
that lie above the flux limit, since the flux limit corre- 
sponds to an intrinsic luminosity that is fainter than that 
for galaxies at higher redshift. If we wish to avoid an 
inherent bias in the pair statistics, we must correct for 
this effect. Furthermore, we must account for any limits 
in absolute magnitude imposed on the sample t o re duce 
the effects of luminosity-dependent clustering (§ |6.l| ). Fi- 
nally, we have demonstrated the importance of specifying 
a limiting absolute magnitude for companions (M2) when 
computing pair statistics. Therefore, we must attempt to 
correct the pair statistics to the values that would have 
been achieved for a volume-limited secondary sample with 
Albright < M < M2. Qualitatively, this correction should 
assign greater importance (or weight) to the more rare 
companions found at the high redshift end of the flux- 
limited sample. To make this correction as rigorous as 
possible, we will use the galaxy LF. By integrating the LF 
over a given range in absolute magnitude, one can obtain 
an estimate of the mean number or luminosity density of 
galaxies in the sample. By performing this integration at 
any given redshift, accounting for the allowed ranges in 
absolute magnitude and the flux limit, it is possible to 
quantify how the mean density varies with redshift within 
the defined sample. This information can be used to re- 
move this unwanted bias from the pair statistics. 



We assign a weight to each galaxy in the secondary sam- 
ple, which renormalizes the sample to the density corre- 
sponding to Mbright ^ M < M2. We must first deter- 
mine A/iini(2i), which gives the limiting absolute magni- 
tude allowed at redshift Zi. At most redshifts, this is im- 
posed by the limiting apparent magnitude m, such that 
Mii^{zi) = m — 51ogdi(z) — 25 — fc(z). At the low redshift 
end of the sample, however, Affaint (defined in § ^.l| ) will 
take over. That is, the limiting absolute magnitude used 
for identifying galaxies in the secondary sample is given 

by 

Mum(z.) = max[Mfaint,m-51ogdL(2:)-25-fc(z)]. (19) 

The selection function, denoted S{z), is defined as the ratio 
of densities in flux-limited versus volume-limited samples. 
This function, given in terms of number density {Sn{z)) 
and luminosity density {Sl{z)), is as follows : 



Ml 

A'/b right 



(j){M)dM 



rM2 



SL{zi) 



J A/brieht 



(j)(M)dM 

ght 

(t>{M)L{M)dM 



(t){M)L{M)dM 



(20) 



(21) 



where L{M) is defined in Equation O. In order to re- 
cover the correct pair statistics, eacn companion must 
be assigned weights WN2{zi) — l/SN{zi) and WL2{zi) = 
1/5*1, (zi). The total number and luminosity of close com- 
panions for the i^^ primary galaxy, computed by sum- 
ming over the j galaxies satisfying the "close compan- 
ion" criteria, is given by = '^jWN2{zj) and = 
J2j ''^L2{zj)Lj respectively. By applying this weighting 
scheme to all galaxies in the secondary sample, we will re- 
trieve pair statistics that correspond to a volume-limited 
sample with Afbright < M < M2. 

6.2.2. Weighting of Primary Sample 

The above weighting scheme ensures that the number 
and luminosity of companions found around each primary 
galaxy is normalized to Mbright < M < M2. However, 
these estimates are obviously better for galaxies at the 
low redshift end of the primary sample, since they will 
have the largest number of observed companions. Recall 
that Nc and Lc are quantities that are averaged over a 
sample of primary galaxies. In order to minimize the 
errors in these statistics, we assign weights to the pri- 
mary galaxies (denoted WNi{zi) and WLiizi)) which are 
inversely proportional to the square of their uncertainty. 
If the observed number and luminosity of companions 
around the i*^ primary galaxy are given by Ni(obs) and 
Li{obs) respectively, and if we assume that the uncertain- 
ties are determined by Poisson counting statistics, then 
= WN2{zi)Ni{obs) ± WN2{zi)^J Ni{obs) and = 

WL2{zi)Li{obs) ± WL2{zi)\/ Li{obs). On average, these 
quantities will be related to expectation values < Nc > 
and < Lc > hy < Nc >= WN2(,Zi)Ni{obs) and < Lc > = 
WL2{zi)Li{obs). Combining these relations yields 



WNl{Zi) 



Ni{obs)wN2{ZtY < Nc> WN2{zi) 



OC WN2{Zt) ^ 

(22) 
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WLliZi) 



Li{obs)WL2{ZiY < Lc> WL2{zi) 



(23) 

That is, the optimal weighting is the reciprocal of the 
weighting scheme used for companions. Therefore, weights 
WNi[zi) = S^izi) and WLi{zi) = SL{zi) should be as- 
signed to primary galaxies. The pair statistics are then 
computed as follows : 



Lr = 



J2i 'WNl{zi)Nci 
J2^'^Nl(z^) 

Y.i^Ll{Zi)Lc, 



(24) 



(25) 



Y.^WLl{Zi) 

It is worth noting that, for a close pair, both galax- 
ies will lie at roughly the same redshift, meaning that 
wi{zi) X W2{zj) ~ 1. We choose not to make this approxi- 
mation, in order to keep these relations valid for pairs that 
are not close, and to allow for future application to pairs 
with additional selection weights. However, we stress that, 
with or without this approximation, the primary weights 
in the denominator provide an overall correction for the 
flux limit, unlike the traditional pair fraction. Note also 
that, for a volume-limited sample, weights for all galaxies 
in the primary and secondary samples are equal, reducing 
these equations to Nc — '^NcjNi and Lc = '^LcjNx, 
as defined in Section h?». 



lies close to the minimum or maximum redshift allowed, 
it is possible that we will miss companions because they 
lie just across this redshift boundary. In order to account 
correctly for this effect, one would need to model the ve- 
locity distribution of companions. As this requires several 
assumptions, we choose instead to exclude all companions 
that lie between a primary galaxy and its nearest redshift 
boundary, provided the boundary lies within /S.v'^^^ of the 
primary galaxy. To account for this exclusion, we assume 
that the velocity distribution is symmetric along the line 
of sight. Thus, as we will miss half of the companions for 
these galaxies, we assign a weight of ti;2„=2 to any com- 
panions found in the direction opposite to the boundary. 
We must also consider how to weight the primary galaxies 
themselves. Clearly, primary galaxies close to the redshift 
boundaries will be expected to have half as many observed 
companions as other primary galaxies. To minimize the 
errors in computing the pair statistics, we assign these 
primary galaxies weights Wy^=0.h. 

To summarize, weights for companions in the secondary 
sample are given by 



Wn-2 = SN^Zi) '^Wh^Wv^, 



(26) 
(27) 



WL2 = SL{zi) ^Wb^Wv^, 

while primary galaxies are assigned weights 

^ SN{zi)wbiWy^, (28) 



6.2.3. Boundary Effects 

A small correction must be made to these weights if a 
primary galaxy lies close to a region of space that is not 
covered by the survey. This will happen if a galaxy lies 
close to the boundaries on the sky, or close to the mini- 
mum or maximum redshift allowed. If this is the case, it 
is possible that close companions will be missed, leading 
to an underestimate of the pair statistics. Therefore, we 
must account for these effects. 

First, we consider galaxies lying close to the survey 
boundaries on the sky, as defined in Section H. For each 
galaxy in the primary sample, we compute the fraction 
of sky with r™'" < r^ < r^'^^ that lies within the survey 
boundaries. This fraction will be denoted /;,. For SSRS2, 
our usual choices of r™'" and r™'*'^ (see § 0) make this 
a very small effect, with /b=l for 99.75% of galaxies in 
the primary sample. Having measured for each galaxy 
in the primary sample, we must incorporate this into the 
measurement of the pair statistics. The first task is to en- 
sure that we correct the number of companions to match 
what would be expected if coverage was complete. We do 
this by assigning each companion a boundary weight wi,2 
= l/Zb, where /h is associated with its host galaxy from 
the primary sample. By multiplying each companion by 
its boundary weight, we will recover the correct number of 
companions. We must also adjust weights for the primary 
galaxies. Following the method described in the previous 
section, we wish to give less weight to galaxies that are 
likely to have fewer observed companions. Therefore, each 
primary galaxy is assigned a boundary weight = ft- 

We now consider galaxies which lie near the survey 
boundaries along the line of sight. If a primary galaxy 



WLi = SL{,Zi)wi,^Wy^. 



(29) 



6.3. Confirmation Using Monte Carlo Simulations 

We will now perform a test to see if this weighting 
scheme achieves the desired effects. To do this, we will 
use flux-limited Monte Carlo simulations, for which the 
intrinsic density and clustering are fixed. Therefore, the 
intrinsic pair statistics do not depend on redshift or lu- 
minosity. If the secondary sample weights are correct, 
the measured pair statistics will be the same everywhere 
(within the measurement errors), regardless of redshift or 
luminosity. We will also check to see if the weights for the 
primary sample are correct. If they are, the errors on the 
pair statistics will be minimized, as desired. 

The flux-limited Monte Carlo simulations were gener- 
ated in a si milar manner to the simulations described in 
Section 5.4; however, a limiting apparent magnitude of 
niB < 15.5 was imposed. Sample sizes of 8000 (South) 
and 4035 (North) were used, providing a good match to 
the overall density in SSRS2. The resulting simulations are 
similar to SSRS2 in all respects, except for the absence 
of clustering. We have already established how the pair 
statistics depend on the choice of Mi and A/2. In the fol- 



lowing analysis, we choose M2=Mi=— 19. In Section S.l 



we outlined reasons for restricting the sample to a fixed 
range in absolute magnitude. Here, we demonstrate how 
the chosen range affects Nc and Lc- For comparison, we 
also compute Nc without normalizing to a specified range 
in absolute magnitude (in this case, wni—wni — ^)- This 
provides some insight into the behaviour of the traditional 
(uncorrected) pair fraction. 
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These tests are most straightforward if the intrinsic pair 
statistics are the same everywhere in the enclosed vol- 
ume. This is not quite true for these simulations, however. 
Galaxies are distributed randomly within the enclosed co- 
moving volume. As a result, the physical density varies 
with redshift as (1 + zf" . In addition, the volume element 
encompassed by the line-of-sight pair criterion Aw varies 
with redshift as (l + z)"^/^ for 50=0. 5. In order to have the 
simulations mimic a sample with universal pair statistics, 
we normalize the sample for these effects by weighting each 
galaxy by (1 + z)~^/^. We stress that this is done only for 
the Monte Carlo simulations. One should not apply either 
of these corrections to real redshift data. 

In Figure |[ the pair statistics are computed for a range 
of Mfaint- In addition, we compute iVc without weight- 
ing by the luminosity function, to demonstrate the dan- 
ger of ignoring this important correction. This statistic 
is directly analagous to the traditional (uncorrected) pair 
fraction used in the literature. It is clear that both 
and Lc are independent of the choice of Mfaint, within the 
errors. This verifies that we have correctly accounted for 
the biases introduced by the apparent magnitude limit. 
In contrast, the unweighted is seen to have a strong 
dependence on Mfaint- As expected, it increases as Mfaint 
becomes fainter, due to the increase in sample density. We 
stress that this does not happen with the normalized Nc 
and Lc statistics, because both arc corrected to a fixed 
range in limiting absolute magnitude. 

Finally, we demonstrate that the weighting scheme used 
for the primary sample (§ |6.2.2 ) does in fact minimize er- 
rors in TVc and L^- Recall that the weighting used was the 
reciprocal of the weights for the secondary sample. Here 
we will ass ume t hat wyi\ cx w-M2{ziY and wi,\ oc w-Li(zi)^ . 
In Section 3.2.2| , justification was given for setting x=—\. 
Here, we will allow x to vary, in order to investigate em- 
pirically which value minimizes the errors. Special cases 
of interest are x=0 (no weighting) and x—\ (same weight- 
ing as secondary sample). The results are given in Fig- 
ure ^ The relative errors in and Lc reach a minimum 
at a; — 1, as expected. Errors are ~ 40% larger if no 
weighting is used (a;=0). For a;=l, errors are much larger, 
increasing by nearly a factor of 5. While errors increase 
dramatically for a; > 0, they change slowly around a;=— 1. 
Clearly, x =—1 is an excellent choice. 

7. A PAIR CLASSIFICATION EXPERIMENT 

The first step in applying these techniques to a real sur- 
vey of galaxies is to decide on a useful close pair definition. 
This involves imposing a maximum projected physical sep- 
aration (r™^'^) and, if possible, a maximum line-of-sight 
rest- frame velocity difference (Aw™'*'^). The limits should 
be chosen so as to extract information on mergers in an 
optimal manner. This involves a compromise between the 
number and merging likelihood of pairs. While one should 
focus on companions which are most likely to be involved 
in mergers, a very stringent pair definition may yield a 
small and statistically insignificant sample. 

In previous close pair studies, the convention has been 
20 hr^ kpc. Pairs with separations of < 
kpc are expected to me rge within 0.5 Gyr (e.g., 
Barnes 19881 iPatton et al. 1997). We note, however, that 



yet to be verified. In earlier work, it has not been pos- 
sible to apply a velocity criterion, since redshift samples 
have been too small to yield useful pair statistics using 
only galaxies with measured redshifts. Instead, all physi- 
cal companions have been u sed, with statistical correction 
for optical contamination (Patton et al. 1997). With a 
complete redshift sample, we can improve on this. This 
can be seen by inspecting a plot of versus Aw for the 
SSRS2 pairs, given in Figure ^. By imposing a velocity 
criterion, we can eliminate optical contamination; further- 
more, we are able to concentrate on the physical pairs with 
the lowest relative velocities, and hence the greatest like- 
lihood of merging. 

We can now use our large sample of \(m-z pairs to shed 
new light on these issues. We will use images of these 
pairs in an attempt to determine how signs of interactions 
are related to pair separation. We begin by finding all 



255 SSRS2 pairs with 



< 



100 /i kpc, computing 



to set r'^'''^ 
20 h-^ 



timescale estimates are approximate in nature, and have 



and Aw for each. Images for these pairs were extracted 
from the Digitized Sky Survey. Interactions were immedi- 
ately apparent in some of these pairs, and the images were 
deemed to be of sufficient quality that a visual classifica- 
tion scheme would be useful. An interaction classification 
parameter (Ic) was devised, where lc=0 indicates that a 
given pair is "definitely not interacting", and Ic=10 indi- 
cates "definitely interacting". In order to avoid a built- 
in bias, the classifier is not given the computed values 
of Tp and Aw. The classifier uses all visible information 
available (tidal tails and bridges, distortions/asymmetries 
in member galaxies, apparent proximity, etc.). Classifica- 
tions were performed by three of us (DRP, ROM, RGC), 
and the median classification was determined for each sys- 
tem. The results are presented in Figure ||. A clickable 
version of this plot, which allows the user to see the corre- 
spondi ng Digitized Sky Survey image for each pair, is avai l- 
able at http://www. astro. utoronto.ca/'^patton/ssrs2/Ic. 

There are several important features in this plot. First, 
there is a clear correlation between 1^ and r^, with closer 
pairs exhibiting stronger signs of interactions. There are 
several interacting pairs with ~ 50 kpc. While 
these separations are fairly large, it is not surprising that 
there would be some early-stage mergers with these separa- 
tions (e.g.. Barton, Bromley, & Geller 1998). An excellent 
example of this phen omenon is the striking tail-bri dge sys- 
tem Arp 295a/b (cf. Hibbard & van Gorkom 1996), which 
has Tp = 95 kpc. However, these systems clearly do 
not dominate; almost all pairs with large separations have 
very low interaction classifications. The majority of pairs 
showing clear signs of interactions/mergers have <J 20 
hr^ kpc . 

There is also a clear connection with Aw. Pairs with 
Aw > 600 km/s do not exhibit signs of interactions, with 
61/63 (97%) classified as Ic < 1. This indicates that in- 
teractions are most likely to be seen in low velocity pairs, 
as expected. We note, however, that there are very few 
optical pairs (i.e., small and large Aw) in this low red- 
shift sample. At higher redshift, increased optical contam- 
ination may lead to difficulties in identifying interacting 
systems when the galaxies are close enough to have over- 
lapping isophotes. Clearly, it is necessary to have redshift 
information for both members of each pair if one is to ex- 
clude these close optical pairs. 
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After close inspection of Figure we decided on close 
pair criteria of r™""" = 20 h^^ kpc and Av"^'^'^ = 500 km/s. 
A mosaic of some of these pairs is given in Figure ^. In this 
regime, 31% (9/29) exhibit convincing evidence for inter- 
actions (Ic > 9), while 69% (20/29) show some indication 
of interactions (Ic > 6). Furthermore, the vast majority 
(9/10) of pairs with clear signs of interactions (Ic > 9) are 
found in this regime. These criteria appear to identify a 
sample of pairs which are likely to be undergoing merg- 
ers; moreover, the resulting sample includes most of the 
systems classified as interacting. 

We also impose an inner boundary of = 5 kpc. 
This limit is chosen so as to avoid the confusion that is 
often present on the smallest scales. In this regime, it is 
often difficult to distinguish between small galaxies and 
sub-galactic units, particularly in merging systems. While 
we are omitting the most likely merger candidates, those 
at separations < 5 h^^ kpc are not expected to account for 
more than ~ 5% of the companions within 20 kpc . 
This expectation, which has yet to be verified, is based 
both on pair counts in i/ST imaging (Burkey et al. 1994) 
and on inward extr apolation of the correlation function 
(Patton et al. 1997). While this inner boundary will lead 
to a slight decrease in Nc and Lc, it should have no sig- 
nificant effect on estimates of merger/accretion rate evo- 
lution, provided the same restriction is applied to compar- 
ison samples at other redshifts. 

8. APPLICATION TO SSRS2 

In the preceding sections, we have outlined techniques 
for measuring pair statistics in a wide variety of samples. 
We have demonstrated a robust method of applying this 
approach to flux-limited samples, accounting for redshift- 
dependent density changes and minimizing differences in 
clustering. We have also selected pair definitions that iden- 
tify the most probable imminent mergers. We will now 
apply these techniques to the SSRS2 survey. As this is a 
complete redshift survey, redshifts are available for all close 
companions; hence, for the first time, we will measure pair 
statistics using only close dynamical pairs. After limiting 
the analysis to a reasonable range in absolute magnitude, 
we compute Nc and Lc for the SSRS2 survey. 

8.1. Defining Survey Parameters 



In Section 3.1 , we emphasized the importance of restrict- 



ing the sample in absolute magnitude, to minimize bias 
due to luminosity-dependent clustering. For SSRS2, we 
first impose a bright limit of Mbright = — 21. All galaxies 
brighter than this are hereafter excluded from the analysis. 
This allows us to avoid the most luminous galaxies, which 
are probably the most susceptible to luminosity-dependent 
clustering; however, this reduces the size of the sample 
by only 0.5%. We also impose a faint absolute magni- 
tude limit of Affaint = —17, which results in the exclusion 
of intrinsically faint galaxies at z <^ 0.01. This guards 
against the possibility that these intrinsically faint galax- 
ies are clustered differently than the bulk of the galaxies 
in the sample. This pruning of the sample is illustrated 
in Figure M. These restrictions allow us to minimize con- 
cerns about luminosity-dependent clustering while retain- 
ing 90% of the sample. The final results are insensitive to 
these particular choices (see Section |8.3|) . 



The above limits in absolute magnitude, along with the 
flux limit, define the usable sample of galaxies. In or- 
der to compute pair statistics, we must also normalize 
the measurements to a given range in absolute magni- 
tude, for both the primary and secondary samples. The 
mean limiting absolute magnitud e of t he primary sam- 
ple, weighted according to Section 6.2.2 , is Mb = —17.8. 
For convenience, we set M2{B) — —18.0 (we will compute 
pair statistics for —19 < M2{B) < —17 in the following 
section). For reference, we note that this corresponds to 
771^ = 15.5 at z=0.017. As we are dealing with a complete 
redshift sample, we set Mi{B)=M2(B) in order to use all 
of the available information. Finally, as we have limited 
the sample using Mbright = —21, this will be used in con- 
junction with M2{B) to derive pair statistics for galaxies 
with -21 <Mb< -18. 

8.2. SSRS2 Pair Statistics 

Using these parameters, we identified all close compan- 
ions in SSRS2. The North sample yielded 27 compan- 
ions, and 53 were found in the South, giving a total of 
80. We emphasize that it is companions that are counted, 
rather than pairs; hence, if both members of a pair fall 
within the primary sample, the pair will usually yield 2 
companions. A histogram of companion absolute magni- 
tudes is given in Figure ^. This plot shows that 90% of 
the companions we observe in our flux-limited sample fall 
in the range —21 < Mb < —18. Hence, galaxies with 
— 18 < Mb < —17 do not dominate the sample. Tables □ 
and ^ give complete lists of close aggregates (pairs and 
triples) for SSRS2 North and South respectively. These 
systems contain all companions used in the computation 
of pair statistics. These tables list system ID, number of 
members, rp {h-^ kpc), Av (km/s), RA (1950.0), DEC 
(1950.0), and recession velocity (km/s). DSS images for 
these systems were given earlier in Figure ^ 

Using this sample of companions, the pair statistics 
were computed. The results are given in Table |l|. Er- 
rors were computed using the Jackknife technique. Re- 
sults from the two subsamples were combined, weighting 
by Jackknife errors, to give Nc = 0.0226 ± 0.0052 and 
Lc = 0.0216 ± 0.0055 x 10^*' /i^L© at z = 0.015. Results 
from the two subsamples agree within the quoted la er- 
rors. 

To facilitate future comparison with other samples, we 
also generate pair statistics spanning the range —19 < 
M2iB) < -17 (see Table |). We note, however, that 
while we account for changes in number and luminosity 
density over this luminosity range (using LF weights de- 
scribed in Section ^ , there is no correction for changes in 
clustering. Hence, our statistics should be considered most 
appropriate for M2{B) — —18, and more approximate in 
nature at brighter and fainter levels. The results in Ta- 
ble H indicate that Nc increases by a factor of 5 between 
M2{B) = -19 and M2{B) = -17, resulting solely from 
an increase in mean number density. The change in Lc is 
less pronounced, with an increase by a factor of 2 over the 
same luminosity range. These substantial changes in both 
statistics emphasize the need to specify M2 when com- 
puting pair statistics and comparing results from different 
samples. In addition, the smaller change in Lc is indica- 
tive of the benefits of using a luminosity statistic such as 
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Lc, which is more likely to con verge as one goes to fainter 
luminosities (see Section 5.4). L^. will always converge 
faster than iVc, thereby reducing the sensitivity to M2. 
Furthermore, it is possible to retrieve most of the relevant 
luminosity information without probing to extremely faint 
levels. For example, for the SSRS2 LF, 70% of the total 
integrated luminosity density is sampled by probing down 
to M2{B) = —18. To first order, the same will be true 
for Lc- Going 2 magnitudes fainter would increase the 
completeness to 95%. While we are currently unable to 
apply pair statistics down to these faint limits, this will be 
pursued when deeper surveys become available. 

8.3. Sensitivity of Results to Survey Parameters 

In this section, we explore the effects of choosing differ- 
ent survey parameters. Earlier in this study, we demon- 
strated that Nc and Lc are insensitive to the choice of 
survey limits in absolute magnitude, provided clustering 
is independent of luminosity and the pair statistics are 
normalized correctly. Here, we test this hypothesis empir- 
ically. 

First, we compute the pair statistics for a range in 
ALfaint, normalizing the statistics to —21 < Mb < —18 
in each case. Figure ^ demonstrates a possible trend 
of decreasing pair statistics with fainter Mfaint- This 
trend, however, is significant only for the brightest galax- 
ies (Mfaint <, —19). This is consistent with the find- 
ings of Willmer et al. (1998), who measure an increase 
in clustering for bright galaxies in SSRS2, on scales of 
Tp > 1 Mpc. For fainter Mfaint, there is no signifi- 
cant dependence. The pair statistics vary by ^ 5% over 
the range —17.5 < Mfaint < —16.5, which is well within 
the error bars. Therefore, we conclude that our choice of 
Mfaint — —17 has a negligible effect on Nc and Lc- This 
implies that, to first order, clustering is independent of 
luminosity within this sample. 

Next, we investigate how the pair statistics depend on 
our particular choices of r^'^^ and Av™'''^, which comprise 
our definition of a close companion. First, we compute 
pair statistics for 10 kpc< r'^'^^ < lOOft,"^ kpc, with 
^^max _ gQQ ]jjjjy'g_ Results are given in Figure |l^. This 
plot indicates a smooth increase in both statistics with 
j,max^ This trend is expected from measurements of the 
galaxy CF. The CF is commonly expressed as a power 
law of the fo rm £,{r, z) = {ro/r)^, with 7=1.8 (Davis & 
Peebles 1983). Integration over this function yields pair 
statistics that vary as r^"'' w r^^, which is in good agree- 
ment with the trend found in Figure From this plot, 
it also appears likely that there are systematic differences 
between the two subsamples. This is hardly surprising, 
since there are known differences in density between the 
subsamples, and it is likely that there are non-negligible 
differences in clustering as well. This cosmic variance is 
not currently measurable on the smaller scales (rp < 20 
kpc) relevant to our main pair statistics. Hence, we 
choose to ignore these differences for now. However, these 
field-to-ficld variations are certain to add some systematic 
error to our quoted pair statistics. 

We also compute pair statistics for a range in Aw™'*'^. 
This is done first for r^^^^ = 20 h^^ kpc, showing the rel- 
ative contributions at different velocities to the main pair 
statistics quoted in this paper. We also compute statis- 



tics using r'^^^ — 100 h^^ kpc, in order to improve the 
statistics. Results are given in Figure 11, Several impor- 
tant conclusions may be drawn from this plot. First, at 
small velocities [Av'"^'^^ <^ 700 km/s), both pair statistics 
increase with Av"'^^^, as expected. This simply indicates 
that one continues to find additional companions as the 
velocity threshold increases. Secondly, it appears that our 
choice of Av"^'^^ was a good one. The < 20 kpc pair 
statistics increase very little beyond Av"^'^^ = 500 km/s, 
while the contamination due to non— merging pairs would 
continue to increase (see Figure^). Moreover, as both pair 
statistics flatten out at around Av"^'^^ = 500 km/s, small 
differences in the velocity distributions of different samples 
should not result in large differences in their pair statis- 
tics. Finally, for r!^^^ < 100 h^^ kpc, the pair statistics 
continue to increase out to ~ 2000 km/s. This indicates 
an increase in velocity dispersion at these larger separa- 
tions. This provides additional confirmation that one is 
less likely to find low— velocity pairs at larger separations, 
thereby implying that mergers should also be less proba- 
ble. 



9. DISCUSSION 

9.1. Comparison with Earlier Estimates of the Local 
Pair Fraction 

All published estimates of the local pair fraction have 
been hindered by small sample sizes and a lack of redshifts. 
In addition, as demonstrated throughout this paper, the 
traditional pair fraction is not a robust statistic, particu- 
larly when applied to flux-limited surveys. The new statis- 
tics introduced in this paper, along with careful accounting 
for selection effects such as the flux limit, yield the first 
secure measures of pair statistics at low redshift. There- 
fore, strictly speaking, the results in this paper cannot be 
compared directly with earlier pair statistics. However, it 
is possible to check for general consistency in results, and 
we will attempt to do so. 

As discussed in Section H, Patton et al. (1997) estimated 
the local pair fraction to be 4.3±0.4%, using the UGC cat- 
alog. The Patton et al. (1997) estimate was based on a 
flux-limited sample with B < 14.5, and a mean redshift of 
z=0.0076. This corresponds roughly to an average limit- 
ing absolute magnitude of Mb = —17.3. Loosely speaking, 
this is analogous to M2. The pair definition used in their 
estimate was < 20 kpc, with no Av criterion. Nc 
may be interpreted as an approximation to the traditional 
pair fraction, provided the relative proportion of triples 
is small. We recompute the SSRS2 pair statistics, using 
^^yinax _ 1000 km/s in an attempt to match the results 
that would be found using no Aw criterion (see Figure pj]) . 
We find Nc = 0.026 ±0.006. This implies a local pair frac- 
tion of 2.6±0.6%. This value is somewhat smaller than the 
earlier result, with larger errors. We strongly emphasize 
that, while these results are broadly consistent, we would 
not expect excellent agreement, due to the improved tech- 
niques used in this study. 

9.2. The Merger Fraction at z ^ 

The quantities Nc and Lc are practical measures of the 
average numbers and luminosities of companions with rel- 
atively high merging probabilities. However, the ambi- 
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guity of redshift space is such that some of these com- 
panions can be entirely safe from ever merging. That is, 
Av = 500 km/s can correspond either to two galaxies at a 
small physical separation with a large infall velocity, or, to 
two galaxies at a separation of 5h~^ Mpc, with no relative 
peculiar velocity. In order to transform from Nc (and Lc) 
to an estimate of the incidence of mergers, we must deter- 
mine what fraction of our close companions have true 3- 
dimensional physical separations of r < 20 kpc. This 
quantity, which we refer to as /sd , has been discussed pre- 
viously by Yee and EUingson (1995). This quantity needs 
to be evaluated based on the small separation clustering 
and kinematics of the galaxy population, tq, 7, and ai2, 
the parameters Av and r™^'^, and the projected separation 
at which two galaxies 'xiptically overlap" in the image. 
The quantity /sd is then evaluated with a triple integral, 
first placing the correlation function into redshift space, 
then integrating over projected and velocity separation. 
For reasonable choices of pair selection parameters, the 
outcome is fairly stable at /3D — 0.5. The most important 
parameter is AV/ai2, which should lie in the range of 2 
to 4. If this parameter is too small, physical pairs will be 
missed; if it is too large, too many distant companions will 
be incorporated. Other reasonable choices include setting 
the ratio of the overlap separation to maximum pair sepa- 
ration to be at least three, and the ratio of the maximum 
pair separation to correlation length to be at least a fac- 
tor of 30. We will take /3d=0.5 to be the best estimate 
currently available. 

We can now estimate the merger fraction (/mg) at 
the present epoch. In this study, we have found 
A''c=0.0226±0.0052. As most companions are found in 
pairs, rather than triplets or higher order N-tuples, this is 
comparable to the fraction of galaxies in close pairs. From 
our estimate of /3D , we infer that half of these galaxies are 
in merging systems, yielding /mg w 0.011. This implies 
that approximately 1.1% of —21 < AIb < —18 galaxies 
are undergoing mergers at the present epoch. We stress 
here that this result applies only to galaxies within the 
specified absolute magnitude limits. Probing to fainter lu- 
minosities would cause /mg to increase substantially. In 
addition, this result applies only to the close companions 
defined in this analysis. Clearly, by modifying this defini- 
tion (and therefore changing the typical merger timescale 
under consideration), the merger fraction would also be 
certain to change. 

9.3. The Merger Timescale 

We now have an idea of how prevalent ongoing merg- 
ers are at the present epoch. In order to relate this result 
to the overall importance of mergers, we must estimate 
the merger timescale (Tmg). We will use the properties 
of our SSRS2 pairs to estimate the mean dynamical fric- 
tion timescale for pairs in our sample. Following Binney 
and Tremaine (1987), we assume circular orbits and a dark 
matter density profile given by p{r) cx r~^. The dynamical 
friction timescale Tftic (in Gyr) is given by 



T, 



fric 



2.64 X lO^r^v. 
MlnA 



(30) 



In A is the Coulomb logarithm. We estimate r and Vc using 
the pairs in Tables || and |[ The mean projected separation 
is Tp 14h~^ kpc. As our procedure already includes a 
correction from projected separation (r^) to 3-dimensional 
separation (r), we take r = rp. Assuming h=0.7, this leads 
to f ~ 20 kpc. The mean line of sight velocity difference 
is Aw ~ 150 km/s. We assume the velocity distribution is 
isotropic, which implies that Vc — V^Av ^ 260 km/s. The 
mean absolute magnitude of companions is Mb ~ — 19 (see 
Figure ||). We assume a representative estimate of the 
galaxy mass-to- light ratio of M/L ^ 5, yielding a mean 
mass of M ~ 3 X 10^°Mq. Finally, Dubinski, Mihos, and 
Hernquist (1999) estimate InA '--^ 2, which fits the orbital 
decay of equal mass mergers seen in simulations. Using 
Equation BOl we find Tfiic ^ 0.5 Gyr. We caution that 
this estimate is an approximation, and is averaged over 
systems with a wide range in merger timescales. Never- 
theless, we will take T^g = 0.5 Gyr as being representative 
of the merger timescale for the pairs in our sample. 

9.4. The Cumulative Effect of Mergers Since z ^ 1 

Now that we have estimated the present epoch merger 
fraction and the merger timescale, we will attempt to as- 
certain what fraction of present galaxies have undergone 
mergers in the past. These galaxies can be classified as 
merger remnants; hence, we will refer to this fraction as 
the remnant fraction (/rem)- We begin by imagining the 
state of affairs at a lookback time oit — Tmg. Suppose the 
merger fraction at the corresponding redshift is given by 
/mg(-z). In the time interval between then and the present, 
a fraction /mg(-2) of galaxies will undergo mergers, yield- 
ing 0.5/mg(^) merger remnants. Therefore, the remnant 
fraction at the present epoch is given by 



Ac 



1 - fmgjz) 
1 - 0.5/mg(z) ■ 



(31) 



Similarly, if we extend this to a lookback time of NT^^g, 
where N is an integer, then the remnant fraction is given 

by 

/rem^l-n /_~n{7^"/^ ' (32) 



0.5/mg(2i) 



where r is the initial physical pair separation in kpc, Vc is 
the circular velocity in km/s, M is the mass (Af©), and 



where Zj corresponds to a lookback time oi t = /Tmg. We 
now make the simple assumption that the merger rate does 
not change with time. In this case, our present epoch es- 
timate of the merger fraction holds at all redshifts, giving 
/mg(-z)=0.011. In order to convert between redshift and 
lookback time, we must specify a cosmological model. We 
assume a Hubble constant of h=0.7. For simplicity, we as- 
sume 90=0.5. Therefore, z = (1 - 3Hot/2)-^/^ - 1. Using 
our merger timescale estimate of T'mg=0.5 Gyr, we can now 
investigate the cumulative effect of mergers. With the cho- 
sen cosmology, z—1 corresponds to a lookback time of ~ 
6 Gyn or 12rmg (A^=12). With this lookback time. Equa- 
tion p3 yields /rem=0.066. This implies that ^ 6.6% of 
galaxies with —21 < Mb < —18 have undergone mergers 
since z 1. 

If the mergers taking place in our sample produce el- 
liptical galaxies, it is worthwhile comparing the remnant 
fraction to the elliptical fraction (cf. Toomre 1977). The 
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elliptical fraction for bright field galaxies is generally found 
to be about 10% (e.g., Dressier 1980, Postman and Geller 
1984). This result is broadly consistent with the remnant 
fraction found in this study. 

While our estimate of the remnant fraction is based on 
our statistically secure measurement of Nc, it also relies 
on fairly crude assumptions regarding the merger fraction 
and merger timescale. In particular, the merger rate has 
been assumed to be constant. There is no physical basis 
for this assumption; in fact, a number of studies have pre- 
dicted a rise in the merger rate with redshift. If this is true, 
we will have underestimated the remnant fraction, and the 
relative importance of mergers. In a future paper (Patton 
et al. 2000), we will address this issue by investigating 
how the merger rate changes with redshift. 

10. CONCLUSIONS 

We have introduced two new pair statistics, N,. and L^, 
which are shown to be related to the galaxy merger and ac- 
cretion rates respectively. Using Monte Carlo simulations, 
these statistics are found to be robust to the redshift- 
dependent density changes inherent in flux-limited sam- 
ples; this represents a very significant improvement over 
all previous estimators. In addition, we provide a clear 
prescription for relating Nc and Lc to the galaxy CF and 
LF, enabling straightforward comparison with measure- 
ments on larger scales. 

These statistics are applied to the SSRS2 survey, pro- 
viding the first statistically sound measurements of pair 
statistics at low redshift. For an effective range in ab- 
solute magnitude of —21 < Mb < —18, we find Nc — 
0.0226 ±0.0052 at 0=0.015, implying that ~ 2.3% of these 
galaxies have companions within a projected physical sep- 
aration of 5 kpc < Tp < 20 h^^ kpc and 500 km/s 
along the line of sight. If this pair statistic remains fixed 
with redshift, simple assumptions imply that ~ 6.6% of 



present day galaxies with —21 < Mb < —18 have under- 
gone mergers since z=l. For our luminosity statistic, we 
find Lc = 0.0216 ±0.0055 x 10^° h^L^. This statistic gives 
the mean luminosity in companions, per galaxy. Both of 
these statistics will serve as local benchmarks in ongoing 
and future studies aimed at detecting redshift evolution in 
the galaxy merger and accretion rates. 

It is our hope that these techniques will be applied to a 
wide range of future redshift surveys. As we have demon- 
strated, one must carefully account for differences in sam- 
pling effects between pairs and field galaxies. This will 
be of increased importance when applying pair statistics 
at higher redshift, as fc-corrections, luminosity evolution, 
band-shifting effects, and spectroscopic completeness have 
to be properly accoimtcd for. The general approach out- 
lined in this paper indicates the steps that must be taken 
to allow for a fair comparison between disparate surveys 
at low and high redshift. These techniques are applicable 
to redshift surveys with varying degrees of completeness, 
and are also adaptable to redshift surveys with additional 
photometric information, such as photometric redshifts, or 
even simply photometric identifications. Finally, this ap- 
proach can be used for detailed studies of both major and 
minor mergers. 

We wish to thank all members of the SSRS2 coUobora- 
tion for their work in compiling the SSRS2 survey, and for 
making these data available in a timely manner. Digitized 
Sky Survey images used in this research were obtained 
from the Canadian Astronomical Data Centre (CADC), 
and are based on photographic data of the National Ge- 
ographic Society Palomar Observatory Sky Survey (NGS- 
POSS). The Digitized Sky Surveys were produced at the 
Space Telescope Science Institute under U.S. Government 
grant NAG W-2166. This work was supported by the 
Natural Sciences and Engineering Research Council of 
Canada, through research grants to R.G.C. and C.J. P.. 
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Table 1 

SSRS2 Pair Statistics 



Sample 


N 


N 

^ ^ comp 


z 




L,(lO^"/i^i0) 


SSRS2 North 


1702 


27 


0.014 


0.0189± 0.0074 


0.0169± 0.0079 


SSRS2 South 


3067 


53 


0.016 


0.0261± 0.0073 


0.0262± 0.0078 


SSRS2 (N+S) 


4769 


80 


0.015 


0.0226± 0.0052 


0.0216± 0.0055 



Table 2 

SSRS2 North : Close Pairs and Triples 



ID 


N 




Aw 


RA(1950.0) 


DEC(1950.0) 


cz 


1 


2 


17.3 


75 


11:18:00.6 


-09:59:44 


1523 


2 


2 


12.7 


385 


12:23:39.7 


-07:24:21 


5183 


3 


2 


19.8 


392 


12:39:42.8 


-05:30:33 


6604 


4 


2 


18.1 


357 


12:49:52.3 


-15:14:45 


4539 


5 


2 


13.0 


183 


12:50:31.4 


-08:55:47 


4149 


6 


3 


11.4 


42 


13:01:21.2 


-11:13:39 


2965 


7 


2 


12.3 


227 


13:17:00.1 


-27:09:24 


1835 


8 


2 


6.5 


198 


13:23:14.6 


-26:24:59 


9596 


9 


2 


18.1 


99 


13:30:31.8 


-15:52:26 


5813 


10 


2 


12.5 


70 


13:32:37.2 


-10:37:55 


6550 


11 


2 


15.7 


298 


13:57:36.5 


-23:03:54 


10592 


12 


2 


12.8 


4 


14:01:52.3 


-24:35:25 


2117 


13 


2 


17.4 


90 


14:10:41.5 


-02:56:41 


1664 
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Table 3 

SSRS2 South : Close Pairs and Triples 



ID 


N 


rp 


Av 


RA(1950.0) 


DEC(1950.0) 


cz 


14 


2 


7.8 


48 


00:13:51.2 


-06:35:57 


8070 


15 


2 


13.7 


276 


00:56:52.5 


-05:04:07 


5771 


16 


2 


18.9 


100 


01:06:26.0 


-28:51:25 


5522 


17 


2 


14.1 


163 


01:10:01.2 


-04:24:30 


5723 


18 


2 


12.4 


160 


01:13:05.7 


-06:51:23 


6365 


19 


2 


12.4 


229 


01:16:22.5 


-17:04:00 


6002 


20 


2 


16.8 


114 


01:34:10.7 


-37:35:09 


5171 


21 


2 


10.7 


173 


02:06:55.3 


-10:22:13 


4143 


22 


2 


6.9 


39 


02:26:44.0 


-11:03:15 


4695 


23 


2 


19.6 


102 


02:55:45.2 


-10:33:21 


4647 


24 


2 


16.0 


78 


03:04:11.5 


-39:13:15 


6188 


25 


2 


19.2 


115 


03:08:51.0 


-09:06:01 


4020 


26 


2 


12.9 


18 


03:22:17.8 


-03:12:17 


2764 


27 


2 


5.4 


43 


04:10:06.5 


-32:59:21 


1085 


28 


3 


7.8 


61 


04:13:09.8 


-28:36:32 


5450 




o 
z 


iU.o 


ZOO 


U4.zy.o4. ( 


-Zv.OL.Zv 


4Uoz 


30 


3 


19.7 


227 


21:08:47.2 


-23:22:03 


10682 


31 


2 


10.8 


99 


21:59:10.7 


-32:13:26 


2674 


32 


2 


11.9 


66 


22:03:22.9 


-28:11:55 


7042 


33 


2 


16.8 


59 


22:19:05.8 


-25:50:09 


6904 


34 


2 


14.2 


107 


22:52:34.7 


-34:09:30 


8700 


35 


2 


19.6 


36 


22:55:33.0 


-04:02:30 


3828 


36 


2 


14.0 


129 


22:57:31.1 


-13:05:22 


3115 


37 


2 


6.5 


47 


23:00:35.1 


-09:15:38 


7360 


38 


2 


19.4 


409 


23:23:42.7 


-39:29:48 


10759 


39 


2 


15.5 


86 


23:45:09.7 


-28:24:56 


8700 



Table 4 

SSRS2 Pair Statistics for Various Choices of M2{B) 



M 


2 








-19 








0082± 





0019 





0135± 





0034 


-18 


9 





0093± 





0021 





0144± 





0037 


-18 


8 





0105± 





0024 





0154± 





0039 


-18 


7 





0117± 





0027 





0163± 





0042 


-18 


6 





0131± 





0030 





0171± 





0044 


-18 


5 





0145± 





0033 





0180± 





0046 


-18 


4 





0159± 





0037 





0188± 





0048 


-18 


3 





0175± 





0040 





0195± 





0050 


-18 


2 





0191± 





0044 





0203± 





0052 


-18 


1 





0208± 





0048 





0210± 





0054 


-18 








0226± 





0052 





0216± 





0055 


-17 


9 





0244± 





0056 





0222± 





0057 


-17 


8 





0263± 





0060 





0228± 





0058 


-17 


7 





0282± 





0065 





0234± 





0060 


-17 


6 





0302± 





0070 





0239± 





0061 


-17 


5 





0323± 





0074 





0244± 





0062 


-17 


4 





0344± 





0079 





0249± 





0064 


-17 


3 





0366± 





0084 





0253± 





0065 


-17 


2 





0388± 





0089 





0257± 





0066 


-17 


1 





0411± 





0095 





0261± 





0067 


-17 








0435± 





0100 





0264± 





0068 
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Fig. 1. — Pair statistics are computed for the volume-limited Monte Carlo simulations. Nc and Lc are given for a range of choices of M2, 
with Mi=M2. Error bars are computed using the Jackknife technique. The dashed lines are the relations predicted using the input SSRS2 
luminosity function, and normalized to match the data at M2=— 16. 
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Fig. 2. — Pair statistics are computed for the flux-limited Monte Carlo simulations. Three pair statistics (unweighted JVc, Nc, and Lc) 3,re 
given, for a range of choices of minimum luminosity Mf^int) with M2=Mi=— 19. Error bars arc computed using the Jackknife technique. The 
unweighted Nc exhibits a systematic increase with Mf^int • This is a consequence of not correcting for the corresponding increase in mean galaxy 
number density. This bias is taken into account when computing Nc and Lc- The horizontal dashed lines match the data at Myuj„j=— 16, and 
demonstrate that JVc and Lc do not depend on Mj^j^j. 
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Fig. 3. — Pair statistics are computed for a range of possible weighting schemes for primary galaxies. Error bars are computed using the 
Jackknife technique. 
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Fig. 4. — Projected separation (rp) is plotted versus rest-fraine velocity difference (Av) for all pairs with rp < 100 kpc. Pairs were 
identified by looking for close companions (real or random) around primary galaxies in the SSRS2 survey. Filled triangles represent real 
companions found in the SSRS2 survey. Open triangles denote random companions found in the Monte Carlo simulations. Note that the 
density of the simulations was matched to the real data sample when generating this plot, ensuring a fair comparison. The upper plot gives a 
wide range in velocity differences, demonstrating how the number of companions in excess of random becomes negligible beyond ~ 1000 km/s. 
The lower plot is a zoomed-in version of the upper plot, giving a better feel for the distribution of physical pairs. 
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Fig. 5.- Pair interaction classification (Ic) is plotted as a function of projected physical separation (rp) and line-of-sight rest-franie velocity 
difference (At)), for 255 unique SSRS2 pairs with rp < 100 h^^ kpc. Symbols are defined as follows : lc=0,l (large open squares), lc=2,3 (small 
open squares), Ic=4— 6 (crosses), Ic=7,8 (small filled triangles), and Ic=9,10 (large filled triangles). High Ic indicates increased likelihood of 
interaction. For plotting convenience, pairs with Av < 10 km/s are assigned Av = 10 km/s. The dashed line marks the close pair criteria 
(rp=20 kpc and Ati=500 km/s) used for computing pair statistics. 
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Fig. 6. — A mosaic of images is given for the 38 close (5 kpc < r-p < 20 kpc) dynamical (Ay < 500 km/s) pairs or triples satisfying 
the criteria used for computing pair statistics. These images were obtained from the Digitized Sky Survey. Each image is 50 kpc on a 
side, corresponding to angular sizes of 1.5'— 10'. 
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Fig. 7. — B— band absolute magnitude is plotted versus redshift for all SSRS2 galaxies with mg < 15.5 and 0.005 < z < 0.05. The curved 
line marks the boundary imposed by this apparent magnitude limit. The upper horizontal line indicates the bright limit imposed (M|-,jigiit = 
—21), while the lower horizontal line denotes the minimum luminosity allowed (Mfaint = Galaxies satisfying all of these criteria (and 

hence used in the calculation of Nc and Lc) are marked with triangles; the remainder are indicated with crosses. 
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Fig. 8. — An absolute magnitude histogram is given for the 80 companions used in the pair statistics. 
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Fic;. 9. Pair statistics are computed for SSRS2 North and South. Nc and Lc arc given, for a range in minimum luminosity Mfajnt, with 
M2=Mi=— 18. Error bars are computed using the Jackknife technique. Both A^c and Lc appear to be independent of Mfg^j^t, within the errors, 
over the range —18 < Mfaint < —16. This implies that, to first order, clustering is independent of luminosity in this regime. 
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Fig. 10. — Pair statistics are computed for Av < 500 km/s, for a range of maximum projected separations (r™*''). A minimum projected 
separation of rp = 5 kpc is applied in each case. Error bars are computed using the Jackknife technique. Both Nc and Lc are cumulative 
statistics; hence, measurements in successive bins are not independent. 
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Fig. 11. — Pair statistics are computed for A?; < 500 km/s, for a range of projected separations. Error bars are computed using the Jackknife 
technique. Both Nc and Lc are cumulative statistics; hence, measurements in successive bins are not independent. 
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